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ABSTRACT 


The air agitated tanks are used in a nunib<‘r of process industries, including metal- 
lurgical industry where these are known as Pachuca tanks. A review of literature 
shows that no generalized correlation is proposed to obtain the mininunn gas flow 
rate as a function of design and processing parameters to obtain complete off-bottom 
particle suspension. In the present work a inatlieniatical model is proposed to predict 
the critical superficial gas velocity (in full-centc'r -column ) with respect to tank as a 
function of important design and process [)<u'<mK'ters like particle diameter, solid con- 
centration, draft tube dimension, etc. Idie validity of tlu' model proposed in prcstmt 
work is checked by comptiriug the results pnrdictt'd by thi.s model with tlu' available 
experimental results obtaint'd by Koich' et al. for bubble; columns with draft tube. On 
comparison it was found that of the t wo approaches ustvl first ap|)roach prt'dicts r(d- 
ativcly better results. .^Iso some modifications in a.pprotich 1 wove doin' and checkt'd 
but the reults did not improve significantly. Tlierc'fon' apprtjtich 1 is t,o be usitd for 
calculating U^c- 
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Chapter 1 


INTRODUCTION 


Gas agitated reactors are used in a number o( process industries, including the met- 
allurgical industry, where they are commonly known as Pachuca or Brown tanks. 
These are widely used for leaching ores in tlie hydrometallurgical production of non- 
ferrous metals like gold, uranium, zinc aiul copper. Another example of the use of 
gas agitated reactors in metallurgical industry is argon stirring of molten steel in 
ladles. However, it is the chemical and hioclu'inical industries where these are used 
extensively. Here these are called bubbU* columns or airlift reactors, and used in 
a wide range of processes such as coal liciui faction, fenmnitatiou, hydrogenation of 
organic compounds, air oxidation ot organic and inorganic compounds, waste water 
treatment, etc. 

Principal objectives in the operation oi gas agitated reactors are the suspension 
of particles, mixing of reagents, mass t.ransler between air bubbles and the slurry. 
The important design parain(d.ers for the Pachuca tanks are flow regiuuj, bubble 
size distribution and coalescence charactt'ristics. a.nd gas holdup. Therefore^ extensive 
studies have been carried out, mostly in laboratory scale bubble columns, to determine 
the effect of different op<?rating parameters such as superficial air velocity, draft tube 
dimensions, tank design and physical prupcrti<'s of the; liciuid/slurry on the process 
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parameters described above. These studies[37] have provided rough guidelines for 
the modification of bubble columns to achieve high process efficiency in different 
applications. However, to date, we still do not have any exact guidelines for the 
design and scale-up of these tanks on a fundamental basis. 

The main objective of gas agitated reactor particle suspension is often very sensi- 
tive to the superficial gas velocity. However, a review of literature indicates that very 
little amount of work has been done on it. No generalised correlation is available to 
predict superficial gas velocity facilitating particle suspension, as a function of var- 
ious design parameters, draft tube dimeusious, particle density, solid concentration, 
etc. In fact, the superficial gas velocity for obtaining complete off-bottom particle 
suspension as a function of design paraiiK'tt.'rs can be predicted. This is done by 
model/correlation developed in this work. 

1.1 Pachuca Tanks 

A stable circulating liquid flow in gas agitat<nl reactors is induced by the apparent 
density difference between the aerated li(iuid (usually in the central zone) in the 
draft tube and the annular region surrounding it. This circulating liquid flow brings 
settling solids in bottom region of the tank into the top region of the lank. Tlu' 
Pachuca tank bottom are usually conical in shape' and a perforated plate i.s used as 
a gas distributor. Conical bottom bubl)h' columns w<'re u.sed by Koide v.i. al.[27]. 
The difference between Pachuca tanks, bubbh' columns and argon stirrc'd ladles lies 
in size and geometry. Pachuca tank diaincirr varices from 5 m to 10 m, diameter of 
argon stirred ladles is usually less then 2 m and that of bubble columns is much less 
than *2 m, usually of the order of 0.15 m. l'h<> lu'ight-to-diameter ratios of Pachuca 
tanks lies between those of argon stirred ladies and bubble columns at values ranging 
from 1.5 to 4. Both Pachuca tanks and bnbbk' columns arc also opcjrated with a draft 



tube positioned in the center just above the perforated plate. Bubble columns with 
draft tube are also known as concentric airlift loop reactors. In Pachuca tanks the 
draft tube diameter to tank diameter ratio is approximately 0.1, whereas in bubble 
columns it is usually greater than 0.4 and at times become as high as 0.87. The draft 
tube serves as an airlift device causing the circulation of slurry in the tank, slurry 
rises in the draft tube and flows downward in surrounding annular space. On the 
basis of draft tube Pachuca tanks can be classified into three kinds, fig.l as discussed 
below: 

1. Full center column type - The draft tube extends from just above the air injec- 
tion point to just below the liquid surface. 

2. Stub column type - In this type draft tube starts just above the air injection 
point and ends just above the conical portion. 

3. Free airlift type - In this type no draft tube is present at all. 

The main advantages in using Pachuca tanks compartjd to mechauiclly agitated 
systems normally used for leaching are : 

1. Less maintenance is necessary due to the absence of moving parts, thus can 
be used for the treatment of corrosive' litpiids and abrasive slurries without the 
problem of corrsion/erosion. 

2. Solids can be handled without any erosion problem. 

3. Aearation of the pulp can be readily provided for enhanced dissolutio:. rates in 
the leaching of gold and uranium. 
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1.2 Research Objective 

In view of our general ability to define the exact nature of two and three phase flows 
without direct observation, it is not surprising to find a similar difficulty in predicting 
pressure drop or energy losses. 

The tentative object of the present work is to construct a model for predicting 
minimum gas flow rate, to obtain complete off-bottom particle suspension in air agi- 
tated tanks valid for a relatively large set of data, parameters being particle diameter, 
particle density, concentration of solids, etc and on this basis to indicate a practical 
approach for design purposes. Full center column air agitated tanks are considered for 
the present work. Then the critical gas velocities obtained that from the e.xperimental 
work of Koide et al.[27] and the model (l<‘ve!op('(l in this work are used to explain the 
effect of various operating and design p<i.rani<?t<>rs such as draft tube diameter, draft 
tube length, particle diameter and concentration of solids. 
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Chapter 2 


Previous Study 


Shah and covvorkcrs have predicted that in inibhle column reactors the hydrodynam- 
ics, transport and mixing properties such <is pressure drop, holdup of various phases 
and interphase mass transfer coefficient de[)cnd strongly on the prevailing flow regime. 
Different criteria have been proposed to differentiate flow regimes. VVallis[49] has 
characterised the upward movement of the bubble swarms into three separate flow 
regimes. These flow regimes occur in ord<u- of increasing gas flow rate. 

Bubbly Flow or Quiescent Bul)bling: 

This regime is characterised l)y bubbl<\s of almost ecjuai size. For this regime 
the superficial gas velocity should be less than .05 m aiul the rise velocity of 
the bubbles in the range .18 to m/s. Tlu'se data should be regarded as a 
guideline for two phase (air an<l waUu') only. 

2. Churn Turbulent Regime or Iletcrogcuius Regime: 

The homogenius gas-in-licpiid dis[)('rsion cannot b(* maintained at higher gas 
flow rates an unsteady flow pattern with chamufiing ocenrs. 'Phis luderogfuiius 
flow regime is characterised by large bubbles moving witli high ris<> velocities in 



the presence of small bubbles. The large bubbles can grow up to a diameter of 
about .15 m. 


3. Slug Flow: 

In small diameter bubble columns, large bubbles are stabilised by the column 
wall leading to the formation of bubble slugs at high gas flow rates. 

Hills[17] extended the flow flow regime maps to high gas and liquid throughputs. The 
dependence of flow regime on column diameter and gas velocity was estimated. Slug 
flow occured in columns of diameter less than .1 m and gas velocity greater than .05 
m/s, churn- turbulent flow occured in columns of diameter greater than .18 m and 
gas velocity greater than .06 m/s. HoniogeniuH flow was possible in all dimensions 
of columns but gas velocity should be low. It is less than .15 m/s for columns of 
dia greater than .15 m and for less it is less than .05 m/s. The type of sparger, 
pysico-chemical properties of liquid, the li([uid velocity can also affect the transition 
between the flow regimes. Like (i)Slug flow can already occur at very low superficial 
gas velocity for highly viscous fluids, (ii) Porous si)argers with mean pore size ies 
than 150 micro m generally produce bubbh' flow up to suijcrficial gas velocities of 
about .05 to .08 m/s. While for pcrforat<td plates with orifice diameter larger than 1 
mm, bubbly flow is possible only at low superficial ga.' velocities. For larger orifice 
diameter bubbly flow may not occur in c<ise of purt> liquids. Thus data above for 
different flow regimes is only approximate. Larger amount of gas can be processed 
through the reactor in bubbly flow. 

In the three phase flow systems, the elhtc.t of solid particles on the flow transition 
has been studied by some investigators (Darton and Harrison, l!J75a,b: Jordani et al, 
1976; Kara, 1981)Darton and Harrison[9, 10] proposeti a trial and error procedure to 
calculate the value of gas and liquid holdup. Where the liquid holdup is empirically 
correlated to gas and liquid velocities and the ratio of clear liquid wake volume to 
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bubble volume. 


2.1 Bubble Dynamics 

Bubble size, bubble rise velocity, bubble size distribution and liquid and bubble ve- 
locity have a direct bearing on the performance of bubble columns. Photographic, 
electrical and optical probe techniques give reliable results only in the homogenious 
or bubbly flow regime. In case the gas is sparged by single orifices, or perforated or 
sintered plates, the diameter of generated bubbles is only slightly dependent on the 
gas velocity. Akita and Yoshida[2] indicates for Dc > 0.3m bubble diameter becomes 
independent of column diameter. Koide et al.[26] reported bubble diameters in a 5.5 
m diameter column and found larger .size.s t.han those predicted by Akita and Yoshida 
and Hughmark[18]. This behaviour was attributed to the relatively large orifice di- 
ameter and insufficiently developed bubbly flow in axial and radial directions. 

Akita and Yoshida[2] found no effect of orifice diameter in bulk region away from 
sparger. Obviously, a balance between coalescnce and breakup rale controlled bubble 
sizes. Also was reported by Deckwer et al.fJl]. But Kumar et al.[28] and Bluivaraju 
et al.[6] proposed correlations for bubbles siz(;s which involve the orifice reynolds no. 
Therefoi'e the effect of the sparger is still not ck'ariflcd. Bubble diameter depciuls on 
the specific gas-liquid system and its proi)erti<\s with respect to coalescnce. 

Calderbank[S] and Vasalos et al.['=10] have reported high rates of bubble coalescnce 
in highly viscous liquids. Bubble brc'akup is due to differences at the interface caused 
by external factors. Low viscosity liejuid .systems are observed to show significant 
bubble breakup. Usually, small bubbh's arc desirable as they yield large interfacial 
ai'cas, which more than offsets decrease in A-/^ with <lecrea.sing diameter. 
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2.2 Gas Holdup 


It can be defined as the percentage by volume of the gas in the two or three phase mix 
in the column. The gas conjunction with the knowledge of mean bubble diameter, 
leads to the determination of interfacial area and in turn mass transfer rates between 
the gas and liquid phase. It depends mainly on the superficial gas velocity, and often 
is very sensitive to the physical properties of the liquid. 

A column diameter greater than 0.15 m is sufficient to obtain the holdup values 
which are close to the ones obtained in larger diameter column, i.e. scaleup has a 
very little effect on the gas holdup. Ueyama and Miyauchi[39] analysis yields that 
the gas holdup in the churn turbulent flow is nearly proportional to Ug and slightly 
decreases with an increase in the column diameter. 

The extreme sensitivity of the holdup to the material system and to the trace 
impurities is not well understood due to large scatter in values obtained by Schugerl 
et al.[36], Bach and Pilhofer[3]. Only a few correlations are based on numerous 
experimental data of Akita and Yoshida[l]. Hikita et al.[l6], Bach and Pilhofer[3], 
and Mersmann[29]. 

Jeffery and Acrivos[21] found that for particles as small as a few ( < 10 ) microns, 
suspensions behave macroscopically as noJl-.N(^wtonian fluids. Therefore, gas holdup 
values involving the non-Newtonian liquids might be useful in analysing the gas-slurry 
systems. Hikita et al.[16] concluded that at low gas velocities the single nozzle gas 
sparger gives lower values gas holdup than a multiuozzle or a porous plate sparger. 

The presence of solids do(?s not alh'ct the gas holdup significantly. Bcgovich 
and Watson[5] combined the available 169 gas holdup points in three phase systems 
and proposed an empirical correlation, ^'ing (h al.[42] applied correlation by Akita 
and Yoshida[l] to t.heir data and conchuled that the correlation given by Akita and 
Yoshida is <iqually adecpiatf; for three jjha.se systems. 



Kato et al.[24, 25], Ying et al.[42], et al.(1980) concluded that an increase in solids 
concentration generally decreases the gas holdup, but the effect becomes insignificant 
at high gas velocities (> 0.1m/ s). 

Freedman and Davidson[13] and Bottou et al.[7] empirically observed insignificant 
effect of the draft tube on the holdup values. On the basis of the other works (Wer- 
land,1978) it can be concluded that the effect of the internals on the gas holdup is 
negligible, and any correlation which fairly represents the values in the presence of 
internals. 

Equations by Akita and Yoshida[l] or llikicta et al.[16] are recommended for less 
viscous and coalescing liquids [viscosity < 0.02 PaS) . These can also be used for 
the three phase systems if the solid phase densities are close to water (~ 1000 — 
1400 Kg/rn^) . 

2.3 Suspension of Particles 

The adequate degree of suspensuion of solid [)a.rticles means hoiuogcuizatioa of slurry 
so formed or a near homogenization such that tlie suspension ensures uniform process 
conditions leading to predictable results. And the adequate degree of mi.xiug is the 
required degree of turbulence which miniinises the mass transfer resistance. 

The suspension is not strictly homogx'uius although appcirently the suspension 
looks complete. The agitation should be such that the solitls arc just suspended and 
the entire solid surface area is e.xpos<'d to tlu' licpud for any reaction in the vessel. .A 
perfect homogenisation at the further expcniditure of iniergy is not desirable. In fact, 
it has been establised that further stirring of an already susp<!aded slurry do('s not 
increases the mass transfer coefficient appreciably. Thus the suspension crit(udon is 
chosen as the condition when no .solid partich's are resting on tlu.- bottom plaXc' of the 
agitated vessel. 
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The small pile of particles remain normally for a very short interval at bottom, and 
at slightly higher speed of liquid, the momentary stay of particles would not exceed a 
second. At this point Narayanan et al.[31] and Zwietering[44] considered suspension to 
be complete. At the point of complete suspension, distribution is fairly homogenious 
when the particles are small and light, whereas coarse and heavy particles remain in 
the lower part of the vessel(Zwietering, 1958). The homogeneity of the suspension is 
not influenced noticeably by the shape of the bottom. 

Regarding homogeneity of the mixture, Durand and Condolios [12] classified par- 
ticles. Particles with a diameter of less than 50 yu m form a homogenious mixture 
even when the movement of the liquid is very slight. Particles from 50 to 200 /i m do 
not form a perfectly homogeneiou mixture, but when completely suspended particles 
were distributed over the whole height of liciuid, it approaches homogeneity at higher 
speeds. Particles above 2 mm move only by successive jumps along the bottom and 
remain in the lower i^art of the liquid. Particles from 0.2 to 2 mm form a transition 
category. 

In systems with low slip velocity, the turbuk'uce generated by this relative velocity 
may be negligible in comparison to the turbulent energy being used up to suspend 
the solid phase and the net effect of the solids is to exert a damping action on the 
fluctuating velocity of the slurry. But an assuini)tion of no slip velocity between the 
particle and the fluid introduces an inherent source of error. 

The only rule, known for scaling up of an agitator which must keep a solid in 
suspension is the rule of e<[ual power [)er unit volume. Zweitering[44] found that 
the required power per unit volume decrc'astjs when increasing the dimensions of the 
apparatus. The rule of equal [)ow('r p<u‘ unit volume is therefore on the safe side. 
Zweitering summed up the conclusion of this section that a confirmation is found 
of the validy of the scaling-up rule based on the power per unit volume, with two 
provisos, one it is conservative because in larg<n' apparatus the power can probabily 
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be lower, and secondly, for radial flow impellers it is more generally valid, in that the 
ratio of tank diameter to impeller diameter can be varied at the same time. 

To effectively transfer matte through a series of stirred tanks, off-bottom suspen- 
sion of the solids in the liquid phase is essential(P.B.Queneau et ah, 1975) An impeller 
imparts both liquid flow and shear stress to a slurry. Flow lifts the solids off the 
tank bottom and disperses, injected or induced gas throughout the tank. Shear stress 
breaks up large gas bubbles and create sufficient turbulence to continually distort the 
bubbles and to prevent a static interface from forming. 

Imafuku et al.[19] observed that the critical gas velocity for complete suspension, 
mostly depends on the liquid near the gas distributor, while Govindrao[14] pointed 
out that the particle diameter and bed volume have a strong influence on the axial 
distribution of solids. 
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Chapter 3 


FORMULATION AND 
APPROACH TO SOLUTION 


In the present study we are trying to construct a model for predicting the critical gas 
flow rate for complete suspension of particles in Pachuca tanks. To obtain this we 
have used the energy balance approach in Pachuca tanks. At steady state the energy 
balance equation in the Pachuca tanks is as given below. 

E, = Ea + E,i+E<ic + Ed. ( 3 . 1 ) 

The terms on the right hand side are energy lo.ss terms and that on left hand side 
energy input terms. E^ is the sole energy source term inside the tank. Ea, is the 
energy loss term in annular region, i.c. region between draft tube and walls of tank. 
Efi is the energy loss in the draft tube. Edc is the energy loss due to sudden contraction 
and the reversal of slurry in lower part of tiuik as the slurry enter? the draft tube from 
the annulus. Edc is the energy lo.ss due to reversal and sudden expansion of slurry, as 
the slurry exits from the draft tube and enters annulus. 

The energy that is entering tlu; tank through flowing gas. is equivalent to the 
amount of work done by gas in expanding isotlnirmally (temperature of slurry and 
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(3.2) 


gas is same) inside the tank, and is estimated by using the following equation. 

E,=m,RT log (h) 

Here rhg is mass flow rate of gas, R is gas constant, T is the temperature of gas at 
which it is flowing into the tank, and Pi is the atmospheric pressure whereas P 2 is 
that at which gas is entering the tank. The assumptions that are implicit in using the 
eq[3.2] are that gas is ideal and whole of the energy lost by gas in expanding under 
isothermal conditions is transferred to the slurry. 

Now there are two different types of flow inside the Pachuca tank, namely, three- 
phase (Gas- Liquid-Solid) flow inside the draft tube and two-phase (Liquid-Solid) flow 
in the annulus. Therefore the energy Ioss('s inside the draft tube and the annulus are 
calculated separately in the following manner. But in both, draft tube and annular 
region flow is similar to closed channel flow. All the friction losses terms used are 
taken from the literature [-16]. The following energy dissipation (per unit mass) term 
is generalised one, for single phase flow in closed channel/pipe. 

P = i(<\/>)2e„ (3.3) 

Here, E is the energy dissipated per unit inass, < V > is average velocity of fluid 
and e„ is the friction loss factor. Now tlu' friction loss factor is calculated in different 
conditions in dififerent way as show Ix'low. 

For sudden contraction. 



e„ = 0.45 

x( l-/i) 

(3.4) 


smaller .1 

:• — section area 

(3.5) 

larger .r 

— seclion area 

for sudden expansion, 


I \ 2 



... = ( 


(3.6) 

and for circular tubes/conduits 

= 


(3.7) 
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(3.8) 


_ .0791 

fle = X V X <1, 


(3.9) 


fl- 

Here, f is friction factor, L is the length of tube/conduit and Rh is mean hydraulic 
radius. Re is Reynolds number, p is density of fluid in the pipe, p its viscosity and 
V its velocity, d, is diameter of pipe in which liquid is flowing. 

The two phcLse flow of water and solid in annulus and draft tube is assumed as 
pseudosingle phase flow. The velocity used in calculating Reynolds number will be 
of slurry. Viscosity and density used in calculations for Reynolds number are that of 
slurry calculated in the following manner. 


pm = p [I - 4>) + Pp <f> (3.10) 

Pm= f) (3.11) 

Here, pm is density of water, Pp is density of solid particles, and pm is viscosity of 
water. The details for the overall density of the suspension given by Eq.[3.10] are 
given in literature[47]. Guth and Simha. [15] had given E(}.[3.11] for concentration 
dependence of viscosity for dilute suspensions of solid spherical particles. 

So as to calculate reynolds no. velocily of slurry is required and since velocity of 
slurry is also not known both slurry velocit.y and reynolds no. are to be calculated 
simultaneously, using the iterative loop met hod. In iterative loop method we assign 
some arbitrary value to one of the two <iuantiti<‘s (in this case velocity of slurry) and 
the other is calculated (Reynolds no.). Now using the calculated value (of Reynolds 
no.) we calculate the other quantity (slurry velocity). Here the loop ends and at this 
point convergence criteria is checked i.<‘. whether it is satified or not If it is satisfied 
we come out of the loop otlujrwise w(> go b<i.ck to the loop. .Now the value of later 
quantity (slurry velocity) is used to calculate <iarlier one (Reynolds no.) and its value 
is used to calculate earlier one (slurry vcl<;city) and so on till the conv<,*rgencc criteria 
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is satisfied, which is checked each time we reaches at the end of the loop. Thus both, 
slurry velocity and Reynolds no. are calculated successively untill and unless the 
values converge in such a manner that successive values of reynolds no. differ by less 
than .01 

We have used the above energy loss/dissipation terms in our approach after mul- 
tiplying these by mass flow rate term, to calculate energy dissipated per unit time. 
The significant energy losses in the annulus are only because of friction between slurry 
and surrounding walls. These are calculated by the following equation. 

Ea = \maX{< V; >)'^ X / ^ (3.12) 

Here, mg is the mass flow rate of slurry in annular region of tank, < V), > is the 
average velocity of slurry in annulus, L is the height of tank, Rh\s radius of lank, and, 
f, friction factor is calculated for tank. The energy losses in draft tube due to friction 
between slurry and walls are calculated using tlie following equation. 

^ rhd X {< K/ >)'^ X / j (3.13) 

Here, rhd is the mass flow rate of slurry in <lraft tube, < Vd > is the average velocity 

of slurry in draft tube, L is the length of draft tube, Rh is the radius of draft tube, 

and f, friction factor is calculated for draft tube. 

At the lower end of the draft tube the slurry enters the draft tube from annular 
region, this leads to sudden contraction of slurry. The energy dissipated/losses due 
to sudden contraction is calculated by th(‘ (‘([nation given below. 

Ed, = ^ (< K: >)■' X [ 0.45 {1-0)] (3.14) 

The above equation is for the case where liciuid is flowing from a larger x-seciion pipe 
to smaller x-section pipe, shown in fig. 2. fl'his equation does not take into account 
reversal losses in Pachuca tanks. Herf‘, < V, > is the slurry velocity at the low(n- end 
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of draft tube, i.e. just inside the draft tube and 

^ area of x — section of draft tube 

/?= - f -- J - ^ 3.15 

area oj x — section oj tank 

At the upper end of the draft tube the slurry comes out of the draft tube and flows 
into the annulus, thereby it expands. The energy losses due to sudden expansion is 
obtained with the help of following equation[]. 


Ede = - rhd X {< Vd >y 


J , 


( 3 . 16 ) 


Here, < Vd > is average velocity of slurry in draft tube. This equation is for the case 
where liquid is flowing from a smaller x-section to larger x-section pipe, shown in 
fig. 2. This equation does not take into account reversal losses in Pachuca tanks. But 
this does not leads to any significant error as is shown later by using the equations 
that take into account reversal losses too. 

On substituting equations 3. 2, 3. 3, 3. 4, 3. 5 and 3.6 in equation 3.1, we get 

[Qu X i< K; '>)'^ X 0.45(1 — /i/) + Qd >< (< 

L 


2Ap,RT CV log ( ^) = \pr. 


1 


+ Qd{<Vd>Y 


L 

Rh 


f] (< K 


Rk 


f 


(3.17) 


Here, Qa is volumetric flow rate in annulus and Qd is volumetric flow rate in draft 
tube, and pm is density of slurry. The above e([uation is rewritten in the following 
form 


t/, 


5C 


2Ap,RT\og (a) 





+ (< v; >fAd [ 0.45(1 -/:()] 


+ f< Ki >fAd 



( 3 . 18 ) 


In final form of energy balance equation < V( >. v('Iocil y of slurry near the bottom 
of draft tube, < Vd >. the average velocity of slurry in draft tube and < Vq >, average 
velocity of slurry in annular region are uuknowu. To calculate these unknowns two 
approaches that have been used by us an' describ(!d in the following lines. 


/ 



'.den expaiision 





I 


-en reduction. 


:1 Coniractioi) in Pipes. 
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3.1 Approach 1 


Rose and Duckworth [33] have developed the most systematically based general cor- 
relation for the prediction of the pressure gradient. Their correlation is applicable to 
vertical, inclined, and horizontal pipe flow of gas-solid and liquid-solid suspentions. 
The correlation developed is based on mechanical energy balance. From their results 
and the basic elements of their correlation, we have taken the minimum transport 
velocity for vertical flows. 

This equation for two phase flows in vertical pipes is used by us because flow 
conditions near lower end of draft tube are very much similar to that in two phase 
flows in vertical pipes. Near the lower end of draft tube flow is nearly two phase row 
as bubbles are not dispersed, and their size is <[uite small. Another reason of using 
this equation is, Imafuku et al. [19] observed that the critical gas velocity for complete 
suspension mostly depends on the liquid flow near the gas distributor. Since draft 
tube and gas distributor, both are centrally positioned in tank draft tube is e.xactly 
over the gas distributor. Therefore the slurry flow near lower end of draft tube ran 
be assumed to be same as near the gas distributor. 

Therefore, Vc, the slurry velocity at the low('r end of draft tube which corresponds 
to minimum transport velocity for off- bottom suspension in vertical pipe flows is 
calculated by following equation 


K = 3.2 X K X li 


0.1 


'Ai 
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XS-^-' X 
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where, K is terminal velocity of particle, Dd is diameter of draft tube, d.p is dian.v'.er 
of particle, R is mass fraction of solids, and s is ratio of solid density to slurry denrny. 
In above equation the terminal velocity of particles is unknown. On equating '.he 
two equations, one of the gravitational force causing the particle to fall down and 
the other of drag force resulting from the fall of particle in fluid, we get an eciua'.ion 



which calculates Vo- 

Vo = 

Cd = — f 1 + 0.1215 ( - 0.05 xiogio (flap) \ 1 ^3 21) 

Rep I- \ p / J 

Equation [3.21] was originally proposed by Beard [5] as a fit to two specific sets of 
data [4, 32] and agrees closely with all reliable experimental and numerical data in 
its range, which is 20 to 260 for particles Reynolds number. A portion of the curve of 
Cd vs NrCp, Reynolds number of particles for spheres is reproduced in Fig. 3. Table 
2 and 3 shows values of Reynolds number and Reynolds number of particles involved 
in our calculations. Also in these tabhjs arc shown corresponding values of Cd, drag 
coefficient. 

Since now we know < >, we can calculate < 14 >? slurry velocity in annulus 

by using the mass continuity equation as shown below. 

TTia = md (3.22) 

pm ^ Qa ~ pm ^ Qd (3.23) 

< v; > x/la = < 14 > xAd (3.24) 

A I 

<Kx> = <14>x-/ (3.25) 

where, rha is mciss flow rate of slurry in annular region of tank and liid is that in draft 
tube. 

Inside the draft tube along with slurry gas i.s also flowing in the form of bubbles. 
As the bubbles rises their size incr{'a.s<‘s and so does slurry velocity in draft tube. 
Therefore, to estimate frictional lossc's in dra.fl tube we calculate < Vd >, average 
slurry velocity inside the draft tube. < 14 > is calculated by the equation given 
below. 

< 14 > = (3.26) 

1 - 9d 
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where, = 1.612 x (3.27) 

(f)ti is the gas holdup in draft tube. In other words, it is defined as volume fraction 
of draft tube occupied by gas. The presence of solids does not affect the gas holdup 
significantly. Begovich and Watson citebib:be combined the available 169 gas holdup 
points in three phase systems and proposed an empirical correlation 

<t,i = (1.612 ± 0.023) X ± ± ° ^ “ (3.28) 

(j)i cannot be calculated directly by eqn[3.26] because Ugc is unknown, therefore 
eqns[3.26] and [3.18] both will be calculated simultaneously in an iterative loop. The 
two are calculated successively till the values converge in such a manner that two 
successive values of superficial gas velocity no more differ than .00001. 


3.2 Approach 2 

It has been observed for baffled tanks, that for small gaps between the turbine and 
the tank bottom {hj H < 1/5) the radial .stream pattern dominates [43]. The particles 
are thrown off towards the periphery, by this radial stream, where they rise along the 
vertical wall of tank. This is equivalent to saying that the solid phase leaves the near 
bottom zone at the points where tlie vertical velocity component of the liquid flowing 
around the particles reaches its maximum value. The value cannot be less than the 
terminal settling velocity of particles. Couscxunuitly it is resonable to assume that the 
minimal flow-around velocity, Vp, near tlu' bottom that is necessary for olf-bottom 
suspension, should be proportional to, Ipi termina.1 settling velocity of particles. 

V;, o k[ (3.29) 

Suspension of slurry by mechanical mixers usually rc‘(|uire a developed turbulent 
regime [Re > 10000) in the tank, (.lonseciinuitly on the tank bottom (or on the 



sediment surface) a turbulent boundary layer is ci'eated with a velocity profile different 
from that in the bulk. At the points where the particles leave the sediment, the 
tangential velocit}'^ component of liquid is substantially less than radial one as the 
measurements indicated [43]. Therefore as a first approximation the stream on this 
path is considered to be one dimensional. Due to the small thickness of boundary 
layer the curvature of the sediment is also neglected. Taking these assumtions into 
consideration Zundelevich et al.[43] expressed the particle flow- around velocity V-p 
near the tank bottom through the relationships defining a one-dimensional flow along 
the flat plate. 

Depending on the particle size, dp, and the viscous sublayer thickness, So, two 
possible versions were cosidered by Zundelevich, namely 


1. If dp < So, the particle is not com])lctcly submerged in the viscous sublayer 
and the average particle flow-around velocity appears to be proportional to the 
particle ratios. 




0.015 X d„ X [' 


(3.30) 


Rcx V 

Here, Vt is the average velocity at the outer border of the turbulent boundary 
layer. 


2. If dp > ^ 0 , the particle is not completely submerged in the viscous sublayer 
beyond which the velocity distribution is characterised by the "’one seventh" 
law. The average particle flow -around velocity becomes 

( (i, X 1/ 

V^ = Ur -oi~ (3.31) 

Here, x is the distance from the flat bottom of th<’ tank. Details are given by 
Zundelevich et ai.citebib:zu. 


Equations [3.29] and [3.30] both includes the unknown value. L't the average 
velocity at the outer border of the turbulent boundary layer, which is close to the 
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mean bulk velocity in the tank U. The ca^e dp > Sq corresponds to the Newton’s 
sedimentation regime. Within the range of action of Stoke’s law and in the almost 
entire transient range dp < 6o . Though in the Pachuca tanks particle size is of the 
order of .074 mm but to match the data available we have taken particle size .198 
mm. Therefore Eq.[3.30] is used in our calculations. We rewrite the equation[3.30] in 
the following form 


Ut 


K + Sit 

(^)' Ref 


(3.32) 


Using the same approximation cis taken by Znidelevich et al. we replace Ut by U in 
equation[3.31] because the average velocity al the outer border of turbulent boundary 
layer is close to the mean bulk velocity, IJ , in the tank. Therefore we get 

Vp + 54 ^ 


U 


(U 




(3.33) 


We further make an approximation, taking Vp, flow-around velocity of particles equal 
to Vof terminal settling velocity of particles. We have used the limiting case since Vp, 
cannot be less than V(>. And this K, is calculated by equation[3.20], earlier used in 
first approach. The final equation that we get is as following 

V; + 54 f- 


U = 


(if)’ rU 


(3.34) 


The above equation calculates minimum, bulk velocity of slurry in the lower part of 
tank for off-bottom suspension of particles. Now we calculate the velocity of slurry 
inside the draft tube using the mass cotitinuily equation [3.24], this we rewrite in 
following manner 


K, = U X (3.35) 

chi 

This is the slurry velocity just entering the draft, tube. To calculate average slurry 
velocity in draft tube we again simultaneously calculate equations [3.18] and [3.27] in 
a loop. Initially some value of is taken arbitrarily and then both are calculated 



successively till the values converge to such an extent that successive values of Ugc 
does not differ by more than 0.00001. 
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Chapter 4 


RESULTS AND DISCUSSION 


Very limited amount of particle suspension data are available in the literature. The 
mathematical models described in the previous section are tested with data from 
Koide et al. [27] work on bubble columns. .All the parameters used in our models 
correspond to those use in his e.xpcriments by Koide et el. In order to ensure the 
validity of the models, results are comparcul lor a large number of parameters such 
as solid concentration particle diameter, draft tube length and draft tube diameter. 
In the figures discussed below in this chapter, the solid or broken lines represent 
the values predicted by mathematical model and symbols represent values obtained 
experimentally (by Koide et al). The solid particles arc of uniform size, with dp = 
0.19S mm and are spherical in shape. 'Flie volumetric ratio or volume fraction of 
solids varied from 0.03 to 0.08 and the size of particles varies in the range 0.1 17 mm 
to 0.498 mm. 


4.1 Solids Concentration 

The increase in concentration of .solids lea<ls to an increase in slurry density and its 
viscosity. This in turn leads to decrease' in slurry velocity. I'herefore to maintain 



minimum slurry velocity for complete olf-bottom particle suspension higher gas ve- 
locity is required. Thus with an increase in concentration of solids the critical gas 
velocity should increase and our models have predicted the same. Fig. 4 shows that 
Ugc values predicted by the two models developed matches quite well with those ob- 
tained empirically. Therefore both the models predict Ugc and its dependence on 
concentration of solids resonably well. 

4.2 Particle Diameter 

Curve 1 and 2 in fig. 5 represents the two different volume fraction of solids .03 
and .08, respectvely. Theoretically[46] the slurry velocity should be just the double 
of terminal velocity of the particles for particle suspension. Empirically it is much 
higher, but has a direct correspondence with Vt, terminal velocity. The slurry velocity 
increases by increasing gas velocity. Vt is directly dependent on the particle diameter 
in such a way that as the particle diameter increases also increases and when 
particle diameter decreases V] also decrease's. Therefore, in other words, Ugc increases 
with increase in particle diameter. Fig. 5 shows that approach 1 has predicted Ugc 
rcsonablly well though the match is not very good. Ugc increases relatively with a 
rapid rate with respect to particle dianu'U'r/tc'rminal velocity. The same figure as 
above shows that approach 2 fails to predict Ugc correctly when the parameter is 
particle cliameter/terminal velocity. 

4.3 Draft Tube Length 

In fig.G there are five distinct sets of data ('acli re|)res<n‘itcd by a different curve, (jurve 
1, 2, 3, and 4 represent different volume fraction of solid in tank. .03, .04, .06, and 
.08, respectively, (.hirve 6 reprc.sent that data in which particle diameter is .117 mm. 



curve 1 approach 1 

curve 1 approach 2 

^ curve 1 expt 
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Effect of solid concentration on Ugc and Comparison of Ugc predicted 
(Lines) by the Two Approaches used and Experimental Data 
(Symbols) from work of Koide et a.l.[27l- 
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log (Terminal velocity) 


Fig. 5 Effect of terminal velocity of particle on U^c ^-ncl Comparison of U^c 

predicted (Lines) by the Two Approaches used and Experimental 
Data (Symbols) from work of Koide el al.[27]. 
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Fig.6 Effect of draft tube length on ( 'gc «iid Comparison of Ugc predicted 
(Lines) by the Two Approaches used and Experimental Data 
(Symbols) from work of Koid'e ot al.[‘i7]. 
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particle density is 4680 Kg/m^ and volume fraction of solids in tank is .04. 

Fig. 6 shows that Ugc first of all deci'eases with increase in draft tube length till 
the length of draft tube becomes equal to the depth of liquid/slurry, (i.e. when the 
upper end of draft tube is in level with the slurry at the upper boundary of tank) and 
after this increases with increase in length of draft tube. The frictional losses in draft 
tube increases with increase in length of draft tube, and this increase in frictional 
losses has to be compensated by increase in energy input rate to maitaincomplete off- 
bottom particle suspension. Therefore higher gas velocity is required. But the values 
obtained empirically shows that Ugc decreases with increase in draft tube length, and 
when draft tube length e.xceeds the liquid level Ugc increases. Till the draft tube 
lies within liquid/slurry the slurry velocity inside the annulus increases, therefore the 
agitation in lower part of tank increases and [^gc decreases. Once it e.xceeds the slurry 
level the velocity does not increase a.iiyniore but the energy losses inside draft tube 
keeps on increasing. Therefore Ugc incnjases to maitaiu particle suspension. Fig. 6 
shows that Ugc values obtained empirically arc ((uite close to those predicted by the 
two approaches. Therefore both the motlels predict the Ugc values as well as their 
dependence on draft tube length rcsonablly well. 


4.4 Draft Tube Diameter 

In fig. 7 there are four different curves 1, 2, 2, and 1 for four different values of volume 
fraction of solid in tank, .02. .04, .06 and .08, resp(’(:tively. Fig.T shows that firstly 
Ugc is decreasing with in(:reas(! in dia.m(‘t<‘r of draft tube and after a certain value of 
diameter of draft tub(? (at which w(' obtain minimum Ug,. ) Ugc starts increasing and 
goes on increasing witli incn’a.se in dianu'tr'r. 

With incrc'ase in diameter of draft tulx^ (k('(!ping tin' tliameL<;r of tank constant), 
the volume fraction of the tank occupicxl by tin* draft tube increases. .Since tluj 
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Fig.7 Effect of diameter of draft tube on U^c and Comparison of Ugc 

predicted (Lines) by the Two Approaches used and Experimental Data 
(Symbols) from work of Koide et al.f’i”]. 
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slurry velocity is expected to be higher within the draft tube compared to annulus, 
increase in diameter of draft tube leads to decrease in Ugc- Another effect of increase 
in diameter of draft tube is increase in slurry velocity in annulus. The gas holdup 
in annular region is mainly determined by the slurry downflow velocity and bubble 
size distribution. If the effective downflow velocity of the slurry in annulus exceeds 
the free rise velocities of the bubbles in stagnant liquid, bubbles are sucked into the 
annular region and swept down by the recirculating liquid. This leads to decrease in 
hydrostatic driving force ( difference in density of slurry in draft tube and annular 
region ) for liquid circulation. Firstly the energy losses due to sudden expansion 
and contraction at the two ends of draft tube decreases with increase in diameter of 
draft tube. When the area of draft tube diameter equals that of annulus the energy 
losses due sudden expansion and contraction reduces to zero. Now when diameter 
is increased beyond this value energy losses at two ends of draft tube due to sudden 
expansion and contraction comes into pla}' and increa.ses with increase in draft tube 
diameter. 

Fig.7 shows that the values predicted by the model developed using both the 
approaches, of Ugc matches quite well witli those obtained ex]5erimentally by Koide et 
al[27]. The values given by approach 1 shows rt^atively better match. The nrinimum 
Ugc is obtained at the ratio of ^ ss 0.0 by Koide et al. experimentally, whereas in the 
case of model developed it is obtained al the ratio of 0.08 for both the approaches. 

4.5 Reversal Losses Coefficients 

The energy loss coefficient e„, for sudden eontracticni and (expansion a.t the draft tube 
bottom and top, respectively, is replacc'd by revc'rsal losses eoelficients at tlie bottom 
and top of draft tube. i.e. AV/ and A'r«, resp<Ttiv(dy. 

AV, = 5.15 X ( ^ ) 


(4.1) 



1 


M 

Aa 


2 



(4.2) 


The correlation for Ktu has been used quite often and also agrees well with exp- 
rimantal data. However, Koide et al. have suggested that ATu depends on the height 
of the liquid level above the draft tube. On the other hand. AV/, depends on the 
configuration near the draft tube bottom. If the geometry of the tank is such that 
it restricts flow from annulus to the draft tube then higher values of Pri (pressure 
loss at bottom of draft tube) should be expected. Therefore a tank with a conical 
bottom (which is the one in consideration), and with a draft tube placed inside the 
cone should yield higher values of Pw compared to a tank with flat bottom (that was 
used by Koide et al. and Miyahara et al.[30]). These were used by Shekhar[3S]. 

The values of Ugc predicted by model using Kru <ind Kri as reversal losses coef- 
ficients are comparable to those predicted without using them. How Ugc varies with 
different parameters and how these valiujs stand against those obtained experimen- 
tally are shown in figs. S to 11. Fig. 12 shows that this model fails in predicting 
how Ugc varies with dralt tube diameter. This is so because Ab; is more sensitive 
to diameter of draft tube compared to A',.,.- Therefore when diameter of draft tube 
increases AT; increases much more than decreases, so energ}’ losses increases. To 
overcome these energy losses more eix'rgy has to be pro\'ided. Thus Ugc increases. 


4.6 Friction Factor Correlation 

In approach one; we had used much more sinq)le but relatively lesa perfect correlation 
for calculating f, friction friction factor. Here we will use a more perfect and also 
more complex correlation for f valid in tlu' full practical range of turbulent flow say, 
3000 < Re < 3,000,000 . 

) = ( 4 xlos( Hi X/-’ )-0.l)'" U.:i) 
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Effect of solid concentration on and Comparison of Ugc predicted 
(Lines) by the Two Approaches used and Experimental Data 
(Symbols) from work of Koide et al.[ 27 ]. 
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Fig.O Effect of terminal velocity of particle on Ugc and Comparison of Ugc 
predicted (Lines) by the Two Api>roac.hes used and Experimental 
Data (Symbols) from work of Koide et aH27]. 
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Fig. 10 Effect of draft tube lengtli on and Comparison of U,jc predicted 
(Lines) by the Two Approaches used and Experimental Data 
(Symbols) from work of Koide et al.[27]. 
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Fig. 11 ElFcct of (liauictcr of draft tulx' on I’^c Comparison of Uyc 

predicted (Lines) by the I'wo .‘\pproachcs used and Experimental Data 
(Symbols) from work of Koidc et al.[27]. 
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The above equation due to Nikuradse[46] provides an excellent fit to essentially 
all reliable data. This equation does not provide us the value of f directly because 
f is also present on right hand side of equation. Therefore it has to be calculated 
inside a loop till two successive values of f do not differ by a constant sufficiently 
small. Kq.[4.3] is based on a form suggested by von Karman(1930) with the constants 
determined from Nikuradse’s data[46]. Figs. 12 to 15 shows Ugc values and how it 
behave.s with various parameters. The Ugc values are nearly the same for both cases, 
i.e. using Eq.[4.3] and Eq.[3.8] for calculating friction factor. But Eq.[4.3] requires 
riKiuire.s lot. more calculations to calculate f, compared to Eq.[3.S] 

4.7 Correlating Single Phase Flow To Multiphase 
Flow 

In t)ur nuMlei we have taken pseudosinglc phase flow. Results have shown that it was 
a. rea.sotuibiy good assumption. Koide et al.[27] have correlated the head losses in 
three phase flow with tlie head losses in single i>hase flow by the following equation. 

113 = f. X hi (4.4) 

< = ( 4 . 5 ) 

ll('r(\ h i is th(‘ head loss in three phase' flow, /q is head loss in single phase flow 
whereas. is volume fraction of licpiid. 

In our model now w(' u.se (h'li.sity and \'is<'osity of water instead of .slurry. .And we 
multiply f. lh(‘ friction 15u-tcu- for singh' pha.sc' by f to (;blaiu friction factor for three 
pha.se flow. I’siiig this we calculate.' all t.lu' values of critical snijerlicial gas velocity 
for otf- bcittom .sus[)ension of particles in Pacliuca tanks. The results arc shown in 
lig.s. Ifi to li). We see iiiat hardly tlu'n' is any significant dilFerenco in the values of 
obt;utie<l tising Koide <*t a), relation and t he one we have assumed. 




Fig. 12 EfFjct of solid roncentratioii on Uyc and Comparison of Ugc predicted 
(Lin<?s) by t,hc Two Ai>proachcs used and Experimental Data 
(Symbols) from work of Koide ct al.[27]. 
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Fig. 13 Effect of terminal velocity of particle on Ugc and Compcirison of 

predicted (Lines) by the Two Approaches used and Experimental 
Data’(SyinboLs) from work of Koide et aL(27|. 
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Fig. 14 


EfF<x:t of <lraft tube length on and Comparison of Ugc predicted 
(Lines) hy the Two Approaciies used and Experimental Data 
(Symbols) from work of Koidc el al.[27]. 




Fig. 15 Eff<‘ct of diameter of draft tube on Ug^ and Comparison of Ugc 

predicted (Lines) by the Two Approaches used and Experimental Data 
(Symbols) from work of Koide et al.[27]. 
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•i of .solid concentration on U^c Compari.son of U,jc predicted 
*.s) by tlie Two Approaches used and Exix'riuu'ntal Data 
hols) from work of Koide ot al.f27]. 
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;.17 Effect of terminal velocity of particle on Ugc and Comparison of V 
predicted (Lines) by the Two Approaches used and Experimental 
Data (Symbols) from work of Koide ot al.[27]. 
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Fig. 18 Effect of draft tube lengtli on Ugc and Comparison of Ugc predicted 
(Lines) by the Two Approaches used and Experimental Data 
(Symbols) from work of Koide ct al.[27]. 
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Chapter 5 


SUMMARY AND COCLUSIONS 


Gas agitat<'cl tanks have been used in tlu* process industry for a very long time. From 
th<‘ engineering point of view the prediction of superficial gas velocity, so as to obtain 
particle suspension in Pachuca tanks, as a function of various design parameters is 
quite important. Still very little htis been done to propose generalised correlations 
predicting tlu' superficial gas vtdocity. Therefore this study proposes a generalised 
correlation prt'dicting suptrrficial gas vcdocity and studies the effect of various op- 
erating and d<;.sign paraincttjrs .such as particle size, solids concentration, and draft 
tube <lim<'nsiuns on sup<‘rficiai gas velocity for off-bottom particle suspension. These 
values of T,;,. predirred by tlu' propos(xl inodel/correlation are compared with those 
obtained I'.xperitne:. tally by Koide et al[27]. 

This stu<iy has tjroposed a generalised and lundaiueutal correlation that helps 
us in under.staiulini: the de[)end<'nce ol superficial gas velocity, facilitating particle 
susp<uiision, on inip'atant chtsign jjaraineters of air agitated tanks. 

1. Increase in concentration of solid particles shouhl lead to decrease in slurry 
velocity near tank bottom. Thendbre particles start settling on the tank bot- 
tom. However, on increasing the superficial gas velocity complete off-bottom 
suspension wiil be restored. 
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2. For a fixed tank diameter, increase in draft tube diameter should lead to increase 
in slurry velocity both in draft tube and annulus. Therefore Ugc decreases, 
minimum Ugc for off-bottom suspension is obtained at ^ = 0.6. Therefore, 
it seems to be somewhat surprising that industrial scale Pachuca tanks are 
operated with ^ ratios around 0.1. 

3. Increase in diameter of particles should lead to increase in terminal settling 
velocity of particle. Therefore slurry velocity near tank bottom is no more suffi- 
cient to keep partich’s suspended. The slurry velocity is increased by increasing 
superficial grus V(docity to a level when off-bottom suspension is attained 

4. Ke«‘ping oilu'r tank dimensions constant, increase in draft tube length leads to 
an increase in slurry velocity. Therefore Ugc decreases, minimum Ugc is obtained 
for c('rtain k'ngtli of draft tube. 

In vi<nv of our gcuu'ral inability to define the exact natue of two and three pnase 
flows withoiit direct t)l).servatiou, we also find a similar difficulty in predicting pressure 
drop accurat<'Iy. From this point of view the results of our models are quite good. 
'I’herefore tin* nuKlels can be used to pix'dict the Ugc to obtain complete off-bottom 
suspension in air agiiateti tanks. 


5.1 Scope of Further Research 

Although th(' correlation proposed is used to predict successfully the empiiical 
velocity ( So as ttj obtain off-bottom particle suspension ) and its dependence tn 
varitius proc<‘ss and <hfsign parameters, tlu? present work prediction could be compa.ed 
with tin* available limited expf'rinumtal data for bubble columns. Therefore it remains 
to be thoroughly checked for Pachuka tanks. Further to establish its comple.ei} 
geiKu-alised form it is essential to compare, tlie predicted Ugc values for large vanaxon 
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on the design and operating parameters used in Pachuka tanks, with the experimental 
data. Because of the difficulties in modelling three phase flows, it is necessary to check 
the reliability of model by conducting experiments on different aspects of particle 
suspension. This e.xercise will also shed some light on the dependence of these process 
parameters on variables that have not been very well understood. 

Therefore the task of determining the superficial gas velocity experimentally for 
complete off- bot tom suspinision as a function of various parameters of Pachuka tank 
and checking the proposed correlation by comparing the values predicted by it with 
thos<' obtained (■\perlm''nl ally is recomended for the future work. 



NOMENCLATURE 


A cr<)ss-.s('rtional area {in") 

, 4,1 cross-sectional area of annular n'gion (in^) 

.4,; cross-sectional area of draft, tube (m"‘) 

Co drag co(‘!Hrient 

D dianu'ter (ni) 

D,j <lianiet<’r of tank (m) 

!),{ (iiauK’ter of draft tube (m) 

(Ij, <iiamet<'r of particles (in) 

I'J euer,e;y flow rate ( lyij;- m/s'’) 

( frii'tion loss factor 
f friction factur 

fy volume fractimi of !i({uid in tank 

g gravitational acceleration (oz/s*) 

h head loss in pii)<‘s (ml 

K energy loss coeilicitmt 

/v',.; enemy loss c<»‘fficient tor re'‘ersal lossi's at, tin* draft tub<“ bottom 

A',.,, eiieryp* loss citefiicient tor reversal losses at, the di'ait tulx* to{) 

A'i, mass t ranfer coeliicimit 

L lens’ th (m| 

dei»th of simrv (m 1 
I., I length of draft tube tin) 

zb m.’is.s flow rate (Kg/s) 

P pri’ssnre ( -N/zzi* 1 

I\ atmosfiheric pressurz* 

P, pn’ssnn* at a [zoint just above (h<‘ gas ditribulor( A/zzz*') 



Pr« r('v<'rsal lost;<\s at tlu' <lraft tube top [lY/in^) 

Q slurry flow rate in tank (;?2^/,s) 

R univ('rsal gas constant (Kg-nj'^/.s'^-K) 

Hh hy<lraulic radius (m) 

T t<>nipcratnr(> (K) 

superficial gas vt'locity hast'd on cross sectional area of tank (m/s) 
b',;,- critical sup<*rficial gas vt'loriiy based on cross sectional area of tank (m/s) 

r bulk slurry vt'iocity in annulus (in/s) 

vt'iuciiy of slurry (lu'sl 
<: \’ ■> averagi' v.'locity of slurry tin/s) 

K,,!) ferininal velocity of parficlt's in stagnant litpud (m/s) 

\ \ vt'locity tif '-lurry in iuinulus (in/s) 

r. \'e!ocily <tf slurry at lliedr.aft tube bottom (in/s) 

\ i; velocity of slurry in drtd’t tidte (m/s) 

(_!reek Symbols 

/> den-dt)' of water 

pm overall dejr-ily of slurry iKg/'tf*) 

/>,, (lensil>' of particles iKg/n/') 

v'» vobmie fracti«»u nf -olids in t;mk 

Oi Viduiiie frai tiou of e.as [diase iii draft tidte 

p vi.'-.<'o;-i!y of Water tKg/iu sl 

/n„ o\'erall vi', Cosily of slurry (Kg/nt-s) 

/' kinematic \i.scosily i s) 

Subsci'ipf s 
.a ammius 

d draft tube 

f liptiid 




slurry 

part 

tank 

top of slurry in tank 
at tliu hn'ol of gas distributor 



1 Tables 



Superficial 

Gas 

Velocity 

(m/s) 

Volume 

Fraction 

Drag 

Coeffi- 

cient 

Particle’s 

Reynolds 

Number 

Reynolds 

Number 

.6990 X 10“* 

.2000 X 10“‘ 

.9666 

116.0 

48040 

.S5S7x 10“2 

.3000 X 10“ * 

.9331 

124.0 

51330 

.99i:ix 10“^ 

.4000X 10“> 

.9108 

129.9 

53800 

.11 lax 10“‘ 


.8938 

. 134.7 

55800 

.12'Jlx 10“‘ 

.6000 X I0“‘ 

.8803 

138.8 

57490 

.i:J26x 10“ ‘ 

.7000 X 10“ ‘ 

.8690 

142.3 

58950 

.M 2 !x l(r‘ 

,8000 X 10“ ‘ 

.8594 

145.5 

60250 

.ir,97x 10“‘ 

.1000 

.8435 

150.9 

62480 

.I7a8x 10“‘ 

.1200 

,8308 

155.4 

64370 

.1908y 10“> 

.1400 

.8022 

159.4 

66010 
























































Ta.ble-2a 


Volume Fraction of Solids = .04 


s. 

No. 

Superficial 

Gas 

Velocity 

(m/s) 

Particle 

Diameter 

Terminal 

Velocity 

Drag 

Coeffi- 

cient 

Particle’s 

Reynolds 

Number 

Reynolds 

Num- 

ber 

1 

.2265 X 10-2 

.1170x10-'' 

— 

.3845x10-' 

1.561 

47.37 

33197.7 

2 

.3509x10-2 

.1360x10-2 

.4496x10"' 

1.327 

63.60 

38350.0 


.6135x10-2 

.1660x10-2 

.5499x10“' 

1.082 

93.20 

46036.3 

4 

.9943x10-2 

.1980x10-2 

.6547x10-' 

.9108 

129.9 

53803.8 

5 

.2050 x 10- ‘ 

.2580x10-2 

.8475x10“' 

.7083 

212.7 

67596.7 

0 

.3068x10-' 

.2980x10-2 

.9743x10“' 

.6190 

277.8 

76428.1 

m 

.5251 x10-' 

.3580x10-2 

.1163 

.5215 

390.0 

89334.6 

8 

.7326x10-' 

.3980x10-2 

.1289 

.4722 

474.6 

97784.7 

9 

.1199 

.4580x10-2 

.1478 

.4135 

616.1 

110311 1 

10 

.1707 

.4980x10-2 

.1604 

.381b 

720.3 

118595.9 1 


50 











Table-2b 


Volume Fraction of Solids = .08 


s. 

No. 

Superficial 

Gas 

Velocity 

(m/s) 

Particle 

Diameter 

Terminal 

Velocity 

Drag 

Coeffi- 

cient 

Particle’s 

Reynolds 

Number 

Reynolds 

Num- 

ber 

n 

.3261x10-2 

.1170x10-^ 

.3974x10-1 

1.461 

53.37 

37403.1 

2 

.5025x10-2 

.1360x10-^ 

.4639x10-1 

1.246 

71.49 

43102.0 

3 

.8759x10-2 

.1660x10-2 

.5666x10-1 

1.019 

104.5 

51620.4 

D 

.1421x10-^ 

.1980x10-2 

.6740x10-1 

.8594 

145.5 

60248.7 

D 

.2962x10-1 

.2580x10-2 

.8718x10-1 

.6693 

237.9 

75615.0 

D 

.4489x10-1 

.2980x10-2 

.1002 

.5850 

310.6 

85480.0 

■ 

.7930x10-1 

.3580x10-2 

.1197 

.4927 

436.3 

99928.8 

8 

.1147 

.3980 ;< 10-2 

.1326 

.4460 

531.0 

109406,5 

S) 

.2118 

.1580x10-2 

.1521 

.3904 

689.7 

123479.1 

LU 

1 .3393 

.1980x10-2 

.1651 

.3603 

800.5 

132800.1 








Table-3a, 


Volume Fraction of Solids = .03 


s. 

No. 

Superficial 

Gas 

Velocity 

(m/s) 

Draft 

Tube 

Length 

(m) 

Drag 

Coeffi- 

cient 

Particle’s 

Reynolds 

Number 

Reynolds 

Number 

1 

.1408x10“^ 

.50 

.8561 

146.6 

60700 

2 

.1086x10-^ 

.85 

.8948 

134.4 

55680 

3 

.1005x10-^ 

1.0 

.9071 

130.9 

54230 

4 

i 

.9223 X 10-2 

1.2 

.9211 

127.1 

52640 


.8587x10-2 

1.4 

.9331 

124.0 

51330 

D 

.1025x10-^ 

1.6 

.9071 

130.9 

54230 

7 

.1031x10-^ 

1.8 

.9071 

130.9 

54230 

8 

.1037x10-1 

2.0 

.9071 

130.9 

54230 

0 

.1041x10-' 

2.1 

.9071 

130.9 

54230 

I iU 

.1044x10-' 

2,2 

.9071 

130.9 

54230 


38 










Table- 3h 


Volume Fraction of Solids = .08 


s. 

No. 

Superficial 

Gas 

Velocity 

(m/s) 

Draft 

Tube 

Length 

(m) 

Drag 

Coeffi- 

cient 

Particle’s 

Reynolds 

Number 

Reynolds 

Number 

1 

.2351 X 10-* 

.50 

.7892 

171.9 

71190 

2 

.1805x 10-^ 

.85 

.8245 

157.7 

65330 

D 

.l66Sx 10-^ 

1.0 

.8357 

153.6 

63630 

4 


1.2 

.8484 

149.2 

61780 

5 


1.4 

.8594 

145.5 

60250 

6 

.1699x 10-^ 

1.6 

.8357 

153.6 

63630 

ra 


1.8 

.8357 

153.6 

63630 

r 

.1721 •< 10“' 

2.0 

.8357 

153.6 

63630 

i 9 

.172GX lU-‘ 

2.1 

.8357 

153.6 

63630 

IBI 


2.2 

.8357 

153.6 

63630 


59 





















Table-4a 


Volume Fraction of Solids = .03 


s. 

No. 

Superficial 

Gas 

Velocity 

(m/s) 

Draft 

Tube 

Diameter 

(m) 

Ratio 

(Ai/A) 

Drag 

Coeffi- 

cient 

Particle's 

Reynolds 

Number 

Reynolds 

Number 

1 

.31S3X 10-^ 

.6200X 10-' 

.4429 

.8325 

154.8 

48474.5 

2 

.2608x 10-' 

.6700x 10-' 

.4786 

.8140 

161.8 

54739.3 

3 

.2088x 10“' 

.7200x 10-' 

.5143 

.7972 

168.5 

61272.8 

n 



.5500 

.7820 

175.0 

68066.7 

5 

.1234X 10-' 

.8200x 10-' 

.5857 

.7679 

181.4 

75113.2 

6 

.9168x 10-'^ 

.8700X 10-' 

.6214 

.7550 

187.5 

82405.5 

7 

.6901 X U)-'^ 

.9200X 10-' 

.6571 

.7430 

193.6 

89937.3 

8 

.5734 X U)-'^ 

.9700X 10“' 

.6929 

.7318 

199.4 

97702.7 

1 

.9'15!) 1U-- 

.1020 

.7286 

.7214 

205.2 

105696.4 

10 

.2477 X H)-' 

.1070 

.76 13 

.7115 

210.8 

113913.4 


60 



















Table-4b 


Volume Fraction of Solids = .08 


s. 

No. 

Superficial 

Gas 

Velocity 

(m/s) 

Draft 

Tube 

Diameter 

(m) 

Ratio 

Drag 

Coeffi- 

cient 

Particle’s 

Reynolds 

Number 

Reynolds 

N umber 

1 

.3562 X 10-^ 

.6200X 10-' 

.4429 

.8188 

159.9 

50066.4 

2 

.2912X 10-^ 

.6700x 10-' 

.4786 

.8007 

167.1 

56535.1 

D 

.2326 X lO-i 

.7200X 10-' 

.5143 

.7843 

174.0 

63281-3 

a 

.iSOOx 10-' 

.7700X 10-' 

.5500 

.7692 

180.8 

70296.2 „ 


.1369x 10"' 

.8200 X 10-' 

.5857 

.7555 

187.3 

77572.0 

6 

.10I5X 10-' 

.8700X 10-' 

.6214 

.7428 

193.7 

85101.6 ; 

B 

.7615 X 10-''* 

.9200 X 10-' 

.6571 

.7310 

199.9 

92878.3 1 

D 

.6307 X U)-'^ 

.9700x 10“' 

.6929 

.7200 

206.0 

100896.3 ‘ 


.1041 X 10-' 

.1020 

.7286 

.7097 

211.9 

109150.0 



.1070 


.7001 

217.7 

117634.3 


Gl 
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